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ABSTRACT 



We present a catalog of photometrically variable stars discovered within two 21.3' x 
21.3' fields centered on the Cygnus OB2 association. There have hitherto been no deep 
optical variability studies of Cyg OB2 despite it being replete with early-type massive 
stars, perhaps due to the high and variable extinction (up to Ay ~ 20) that permeates 
much of the region. Here we provide results of the first variability study with this 
combination of spatial coverage (~ 0.5 deg) and photometric depth (R ~ 21 mag). We 
find 121 stars to be variable in both R- and /-band, 116 of them newly discovered. Of the 
121 variables, we identify 27 eclipsing binaries (EBs) and eclipsing binary candidates, 20 
potential Herbig Ae/Be stars, and 52 pulsating variables. Confirming both the status 
and the cluster membership of the Herbig Ae/Be stars would address the uncertainty 
regarding the age and star formation history of Cyg OB2. We match our catalog to 
known variables and binaries in the region, 2MASS near-IR (NIR) data, and Chandra 
X-ray observations to find counterparts to new variables in other wavelengths. 

Subject headings: Cyg OB2 — binaries: eclipsing — stars: massive — stars: variables: 
general 



1. Introduction 

Procuring accurate fundamental measurements of the masses and radii of stars remains one 
of the most outstanding questions in astronomy. This is crucial in the high-mass regime, which 
is essential for determining the slope of the initial mass function (IMF) for high-mass stars. The 



fundamental determinations of stellar masses with errors <1— ! 

reject models with incorrect underlying physics. As emphasized in the review by 



Torres et al. 


2010) in order to 


in the review by 


Torres et al. 



1 Based on observations obtained using the 1.3m McGraw-Hill telescope at the MDM Observatory 

department of Astronomy, The Ohio State University, 140 W. 18th Ave., Columbus, OH 43210, USA 

3 Hubble and Carnegie-Princeton Fellow at Carnegie Observatories, 813 Santa Barbara St., Pasadena, CA 91101, 
USA 
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(|2010l ). there are only 95 detached binary systems of any mass in which the mass and radius of 
both components are known to ±3% or better. All but one are eclipsing binaries, with the sole 
exception being the Sun's nearest neighbor, a Centauri. Of the 95 total systems there are only 
eight with M > 15 M , three with M > 20 M , and none with M > 30 M Q . Even when applying 
a less stringent uncertainty cutoff, only 14 stars with M > 30 M & have both their mass and radius 
known to ±10% or better (jBonanosI 120091 ) . This paucity of accurately determined stellar masses 
and radii at the high-mass end requires an expedient redressing. Fundamental measurements of 
properties for a statistical sample of massive stars will allow us to constrain theoretical models, 
anchor the IMF at the high-mass end, and better comprehend stellar structure and evolution. 

The Cygnus OB2 association presents the ideal laboratory within which to learn about the 
fundamental properties of massive stars. It is the nearest massive star-forming region, making it 
a prime target observationally, even in spite of the extinction in the larger Cygnus region, which 
can range up to Ay ~ 20. Furthermore, it is known to be replete with an estimated 2,600 ± 400 
OB stars, includin g 120 ± 20 O-type stars, and have a total inferred mass of - (4 - 10) x 10 4 M 



( Knodlseder 



200 



(j). The region has bee n the focus of several studies at a variety of wavelengths, 



dating back to Munch MorgarJ ( 19531 ). who initially noticed a clustering of 11 blue giants near 7 
Cygni. Ijohnson Sz Morgan! ( 1954! ) were the first to obtain accurate spectral types, finding 10 to be 
O-type and one to be an early B-type star. The subsequent discovery of the heavily reddened and 



extrem ely luminous star Cyg OB #12 ([Morgan et al.lll954af) w as followed by seven ([Morgan et al 



1954b! ) and then 12 additional early- type stars (|Schultdll956bl ). The first comprehensive study of 
the morphology of the association estimated more than 3,000 total me mbers, of which at le ast 300 
are early-type (OB) stars, and a total mass of ~ (0.6 — 2.7) x 10 4 M & ( Reddish et al.lll966l ). 



The first synoptic optical study of VI Cygnus produced UBV photometry of 799 stars in the 
region and obtaine d spectral types for 76 stars , most of which are the brightest and bluest members 

They found a distance of 1.7 kpc and the slope of 



1991 



of the association (jMassev Sz Thompson 
the IMF to be T = —1.0 db 0.1. This provided the first photometric catal og of Cyg OB2 and spec- 
troscopically identified the largest number of OB association members. iPigulski fc Kolaczkowski 
(|l998l ) performed the first variability survey of the region, obtaining /-band photometry of the 
central 6 x 1 2 arcmin 2 . They fo und 29 variables, five of which are O-type, and at least five eclips- 
ing binaries. Knodlseder (|200ol ) studied the region using 2MASS data to characterize the stellar 
population. They aver an IMF slope of F = — 1.6 ± 0.1 and assert that the large total mass, 
high central density, and overall compactne ss of the region requi re that it be classified as a young 
globular cluster, a claim originally made bv lReddish et al.l ([1966). 



Precise classification notwithstandin g , early -type members continued to be discovered spectro- 
scopically in the region. IComeron et al.l (120021 ) used 2MASS photometry to identify hundreds of 
potential early-type stars and spectroscopically found 77 of them to be early-type members. Of 
those, 31 had prior classifications, 24 of which are O-type and the remaining seven B-type, causing 



them to extrapolate that a total of 60 of t heir sam 



contains between 90-100 total O-type stars. lHanson 



p ie of 77 are O-type stars and that Cyg OB2 
(|2003l ) took spectra of 14 of the remaining 46 
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and identified all 14 to be early-type, perhaps indicati ng that t he rem aining 31 (one was discovered 
to have been previously studied) are also early-type. iHansonl (j2003l ) then extrapolated the locally 
derived cluster luminosity function to argue that ther e could be tens or ey en one hundred similar 
massive open cl usters within the Milky Way . Recently, iKiminki et al.l (|2007l ) built off the canonical 
optical study of lMassev Thompson! (|199ll ). spectroscopically following up 183 of their photomet- 
rically observed Cyg OB2 members. Seventy-three possessed extant spectral classifications as OB 
stars and 73 of the remaining 110 were also found to be OB stars. Among these they flagged 36 
probable and nine possible single-lined spectroscopic binaries (SBls) as well as three new and eight 
candidate double- lined spectroscopic binaries (SB2s). 

Several addition al binaries have been discovered in the Cyg OB2 region in the past few years. 
Kiminki et al.1 (|2008l ) presented orbital solutions for two SBls and one SB2 that they had identified 
as SB-candi dates in their earl ier work. Cygnus OB2 #9, a non-thermal radio-emitter, was confirmed 
as a binary (|Naze et al.ll2008h . supporting the theoretical wind- wind collision s cenario. Two EBs i n 
the Cyg OB2 r egion were found via data in the NSVS survey (|oteroll2008«] fl Istroud et allhoioh . 
Kiminki et al.1 (120091 ) determined the orbital parameters of four additional OB SBs. This brings 
the total number of known OB binary systems with full photometric and/or spectroscopic solutions 
within Cyg OB2 to 18. 

Furthermore, the uncertainty regarding the age and star-formation history of Cyg OB2 makes 
this cluster all the more int riguing. The abundance of heavily-reddened early-type stars found by 
Massey &: Thompsonl ()199ll ) initially suggested a single epoch of star fo rmation and a c o rrespo nding 
age of 1-3 Myr. Several of the e arly-t ype candidates identified by IComeron et al.l ([2002]) were 
confirmed as OB stars by IHansonl (|2003l ). but their dubious sta tus as cluster m embers casts doubt 
on the 1—1.5 Myr age that their membership would necessitate. IHansonl (12003) goes on to question 
the membership status of the remaining stars in the I Comer on et al.l (|2002l ) sample, and the presence 
of evolved O-stars such as the s upergiant Cyg OB2 #9 requires an age of at least 3 Myr. More 
recently, IPHAS observations by iDrew et al.l (|2008l ) revealed a clustering of ^200 A-type dwarf s 
primarily to the south o f Cyg OB2 and withi n the 5-7 Myr age range. I Comer on et al.l (J2008) 
used the same method of IComeron et al.l (|2002j ) to investigate an annular region between 1 and 2 
degrees whose center lies slightly to the north of Cyg OB2. They uncovered a diffuse spread of 
evolved early-t ype stars wtih ages up to ~10 Myr that they identify as likely being field population 
stars. Finally, I Wright et al.l (|2010l ) used isochrone fits to NIR color-magnitude diagrams (CMDs) 
to derive ages of 3.5 and 5.25 Myr for their two fields. Given the large spread in ages, it is unclear 
whether the central core of Cyg OB2 is merely the latest in a series of epochs of star formations, 
if it is the product of a parent cluster, or if the populations of older and more evolved stars are 
simply not cluster members. 

The Cygnus OB2 association is also interesting as a source of high-energy (>MeV) emission. 
Due to limitations of angular resolution, most observations thus far cannot clearly separate Cyg 
OB2 from nearby parts of the Cygnus region more generally. This region is bright in gamma 
rays from both point sources and diffuse emission, ranging from the MeV to the GeV to the 
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TeV (jHartman et al.1 Il999l ; IProdanovic et al.1 120071 : lAbdo et al.1 120071 ; lAlbert et all |200J; iRomero 
2010) . Neutrino s and gamma rays from this region will be important for testing models of cosmic 
ray production ( Anchordoqui et al.l 120071 ; iBeacom &; Kistlerl 120071 : lHalzen 6 Murchadhal 120071 ; 
Anchordoqui et al.ll2009i ). With present and future instruments, it should be possible to better 
resolve the young Cyg OB2, which would allow the isolation of the high-energy emission from 
pre-explosion massive stars as opposed to post-explosion remnants. 

This paper describes the findings of the first deep optical two-band (Rc, Ic) photometric 
variability survey of the Cygnus OB2 region. We find 121 variables and provide variable classi- 
fications for them. In §2 we give a full description of the photometric observations. The details 
of the photometric calibration and light curve production are in §3. We present the results of 
the variability search in §4 and discuss the different variability classes separately. We match our 
variability catalog to the 2MASS Point Source Catalog in ord er to find NIR counterparts and then 
to a recent X-ray survey of Cyg OB2 (jWright &; Drake! l2009i ) . both of which are discussed in §5. 
Finally, §6 summarizes the results of the survey and future directions of the project. 



Data 



2.1. Observations 



The data were obtained on the MDM 1.3m telescope using the 4K imager, which has a field 
of view (FoV) of 21.3' x 21.3' and a pixel scale of 0.315" per pixel. The observations consist of 60s 
exposures taken through the Cousins i?-filter and 10s exposures taken through the Cousins /-filter. 
The data were collected over 19 nights during the fall of 2007, spanning September 12 - October 
24, and 18 nights during the summer of 2009, spanning May 30 - June 25. Two pointings were 
used to provide greater spatial coverage of the Cygnus OB2 region. The first field was centered 
at (a = 20 h 33 m 17 s , 5 = +41 d 14 m 15 s ) and the second field to the north at (a = 20 h 31 m 53 s , 5 = 
+41 d 27 m 00 s ), both chosen in order to cover the central core of the cluster and coincide with the 



study of lMassev &; Thompson! ()199ll ). TableQ]gives an overview of the observations. Figure Q] shows 
a 2 x 2 deg 2 image of the Cygnus OB2 region with overlays of th is study, the prev i ous v ariability 
study of lPigulski fc Kolaczkowskil (j 19981 ). and the X-ray survey of lWright fc Drake! (J2009J). 



2.2. Basic Data Reduction 

We apply basic CCD reductions using a custom software pipelin^] developed by Jason Eastman. 
The routine performs overscan subtraction, trims the images, and uses dome flats to flat field the 
data. The 4K chip suffers from crosstalk contamination, a subtle but insidious effect, and so the 



1 http://www.astronomy.ohio-state.edu/~jdeast/4k/proc4k.pro 
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procedure was modified to properly correct for this. Crosstalk results from the chip being divided 
into four quadrants of equal pixel size that are read out by four different amplifiers simultaneously. 
During the readout phase, signal from bright stars being read out on a given quadrant, referred 
to as the source quadrant, will reverberate through the electronics and leave residual signal on the 
other three quadrants, called victim quadrants. The source signal is reflected across the three axes 
of symmetry of the chip, making the contamination predictable and predominantly tractable. 

Before the residual signal can be removed, a set of crosstalk coefficients, which specify the 
fraction of the source signal that is reflected onto each victim quadrant, must be calculated. The 
coefficients are small enough (~ 0.04% — 0.1%) that the problem can be sufficiently dealt with 
via a simple solution. Using a set of a few hundred images of our data from both the 2007 and 
2009 seasons, we find that the coefficients are not only constant in the count regime 45000 < 
counts < 65535, but also stable across the two-year baseline. Information cannot be backed out of 
reflected signals due to saturated (65535 counts on the 4K imager) sources, so pixels contaminated 
by saturated stars are set to zero, removed from the light curves, and otherwise ignored. 



3. Photometric Calibration and Light Curve Production 



packa ge (Stetson 



ISIS package (Al ard &; Lupton 



(2003) and 



1987, 



Hartman et al. 



199JL the procedural details of which are discussed in iMocheiska et al 



1998 



Alard 



We use as trometry.net ( Lang et al.ll201C ) to d o astrometry and we extract photometry using the 



2000) in conjunction with the DAOP HQT/ALLSTAR 



(|20Q4j ) . We did not observe standard stars, so this requires that we 



match to other calibrated catalogs and apply a computed magnitude offset. For the i?-b and data we 



matc h our full catalog to an overlapping region of stars observed by the IPHAS survey (IDrew et al 



20051 ). We uniquely match six stars from their catalog to ours using a matching radius of 1". We 
then apply the transformation equation given by Lupton et al. (2005)0, 



R 



0.2936 • (r 



0.1439, 



(1) 



to transform the IPHAS SDSS magnitudes into Cousins i?-band magnitudes. From this we derive 
th e magnitude offset and ca li brate our i?-band data. To calibrate our /-band data with the catalog 
of lPigulski &: Kolaczkowskil (119981 ) we use a more stringent star selection, as there is more spatial 
overlap between the two data sets. For a star to be used for the calibration, it has to satisfy the 
following criteria: 1) have a unique match within l", 2) have no additional matches within 5", 3) 
not be classified by them as a variable, and 4) not be saturated or blended in our own database. 
This ultimately leaves 133 stars that we then use to calibrate our /-band data. As each of our 
fields has a different reference frame, there is a slight magnitude offset between fields that can be 
calculated using stars in the overlap region. 



2 http://www.sdss.org/dr6/algorithms/sdssUBVRITransform.html 



- 6 - 



We perform a variability search using the VARTOQL S light curve analysis program (Hart man et al 



2008). We employ both Analysis of Variation (AoV) flS chwarze nberg-Czernv 



and Lomb-Scargle (LS) (|Lomblll976l ; IScargldll982l ; LPress fc Rvbicki 



198 



9; 



2005) 



1989; Devor 
Press et al]ll992h period 



search algorithms. We then inspect light c urves with high signific ance by eye to find periodic vari- 
ables and use the Stetson variability index ([Welch Stetsonlll993l ) to select irregular variable stars. 
Figure [2] shows the light curve rms as a function of reference magnitude for our 23,561 i?-band light 
curves and 12,611 /-band light curves. We find variables as deep as R ~ 21 and I ~ 19.5 and go 
fainter in R due to the longer exposure time. Many stars with high rms are not listed as variable 
because we establish the criterion that stars must exhibit simultaneous two-band variability, which 
eliminates false positives and mitigates a ny pote n tial r esidual effects of crosstalk contamination 
from saturated sources. As described in IStetsonl (|1996l ). stars displaying variability in only one 
band will have a variability index near zero. Only stars that exhibited clear two-band variability 
are selected, 121 in total. 



Variable Stars 



Prior to discussing the parame ters of the variables, the various naming conventions must be 



expl ained. 
that 



Morgan 



et al 



and ISchultel |l956b 



Mun ch fe Morgan! ()1953l ) numbered the first 11 stars discovered in Cy g OB2, a system 



1954al ) expanded for star number 12, iMorgan et al.l (ll954bl ) for stars 13-19, 



|aj) for stars 20-85. We follow the convention of the literature, which is to 
precede these stars with 'Schulte'. All variable stars in t he G CVS are identified as such. The 



primary numbering system is that of lMassev Thompson! (| 199 11 ), for which the designation MT91 
precedes star ID numbers, from 1-799. The variability study of iPigulski &: Kolaczkowskil (|1998l ) 
extends this system to 975, and is referred to as MTE (MT91 Extended). The designations A36 
and A45 (GS C 03161-00815, 2 MAS S J20294666+410 5083, respectiv e ly) ar ise from the convention 
introduced by IComeron et al.1 (|2002j) and followed by lKiminki et al.l (|2009i ). Stars within our own 
variability catalog are given a VarlD, and numbers increase from brightest to faintest mean .R-band 
magnitude. 



Pigulski Kolaczkowskil (|1998l ) performed the first variability survey of the region, detecting 



455 stars in /-band and finding 29 of them to be variable, including five EBs and five O-type stars. 
Our study goes much deeper, covers significantly more area, and furthermore provides a two-band 
photometric catalog. We identify 121 stars that exhibit variability in both R- and /-band. Of 
these, 10 are EBs, nine of them newly discovered, 17 are EB-candidates, 20 are potential Herbig 
Ae/Be stars, and 52 are pulsating stars. We first discuss the results of matching our catalog to the 
literature. Then we discuss each variability class individually: first the EBs and EB-candidates, 
then the pulsating variables, and finally the potential Herbig Ae/Be stars. 
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4.1. Literature Comparison 

We find 121 stars to be variable in both R- and /-band within the Cygnus OB2 association. 
After doing a literature search, we match our catalog to all previously known variable stars in Cyg 
OB2, listed in Tabled using a matching radius of 1". Of the 30 known variable stars in Cyg OB2, 
we recover variability in at least one band for nine of them. Seven stars are without matches while 
the remaining 14 are either saturated, photometrically constant, or outside our FoV. We discuss 
those with matches here. 



V2186 Cyg is an EB of the Algol type with a period of 2.9788 days (jPigulski fc Kolaczkowski 



19981 ). Though it is saturated in our it!-band data, we match to it in /-band and successfully 



recover the period. MT91 460 is an irregular variable with an /-band amplitude of 0.28 mag 



(jPigulski Sz Kolaczkowskilll998l ). We do not match to it in our i?-band data, but find it to exhibit 
a 0.7 mag drop in /-band, with the overall shape of the light curve resembling that of a long period 
variable. Schulte 57, VarlD 020 in our catalog, was previously classified as a variable of unknown 
type with a 0.728 day period. We recover the period in both bands and classify it as a pulsating 
star. 



The remainder of the stars discussed here are part of the MTE system, introduced bv lPigulski Kolaczkowski 



(|1998l ) to extend the catalog from 799 to 975 stars. They report MTE 831 to be irregular with an 
/-band amplitude of 0.22 mag. It matches to our VarlD 050, with a 10 day period and a 0.3 mag 
amplitude in both R- and /-band. MTE 849 (VarlD 021) is an EB of the Algol type. We report a 
period of 2.89 days, the first known for this system. MTE 876 (VarlD 036), previously a variable 
of unknown type, appears to have a characteristically pulsating light curve, with a period of 1.53 
days and a 0.1 mag amplitude. We agree with their classification of MTE 900 (VarlD 038), which 
varies erratically with an amplitude of 1.0 mag on long timescales, perhaps indicating that it is 
an irregular or long-period variable. MTE 916 (Var lD 011) similarly appear s to be a long-period 



variable, though with a smaller 0.3 mag amplitude. iPigulski &: Kolaczkowskil (|1998l ) classify MTE 
973 as a potential EB, and we find it to be constant in /?-band but to have a 0.2 amplitude and a 
0.88 day period in /-band, though with no obvious eclipsing signature. 



4.2. Binary Systems 

In addition to the eight previously known EBs within Cyg OB2, there are 13 OB SBs in Cyg 
OB2, four of which are single-lined and 11 of which are double-lined. It should be noted that a 
few of these spectroscopic systems are also definite or possible eclipsing systems. Those that are 
confirmed EBs, Schulte 27 and A36, are treated here as photometric variables. Those that are 
unconfirmed, MT372 and Schulte 3, are included here only as SBs. Of the 121 stars we identify 
as photometrically variable, 10 are EBs and 17 are EB-candidates. Figure [3] shows sample light 
curves for four of the EBs and figure 2] shows four of the EB-candidates. 
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To determine if any of our variable stars are SBs, we mat ch to all confirmed an d potential SBs 
found in the literature, most of which are from the catalog of lKiminki et al.l (120071). They present 
spectr oscopic data for 146 OB stars within Cyg OB2, all initially discovered bv lMassev fc Thompson 
(|199ll ). and find 73 new OB stars. Restricting themselves to the 120 stars with the best data, they 
identify 36 probable and nine possible SBls and three new and eight possib le SB2s. Table p lists 
confirmed and potential SBs within Cyg OB2. Given that the catalog of iKiminki et al.1 (|2007l ) 
reaches a limiting magnitude of V ~ 15 and our upper limit is R ~ 13.5, the vast majority of the 
stars are saturated in our catalog, and thus we find no matches. 



4.3. Pulsating Variables 

We find 52 of our stars to be pulsating variables, with i?-band amplitudes ranging from 0.1 to 
0.6 magnitudes. Figure [5] shows sample light curves of four of the pulsating variables. The periods 
range from ~0.05 to ~60 days. A histogram of the period distribution from to 20 days is shown 
in figure [H 



4.4. Herbig Ae/Be Stars 



Initially identified by iHerbia (|1960i ). Herbig Ae/Be stars are pre-main sequence (PMS) stars 
with masses > 2 Mq and are the intermediate-mass analogs of T Tauri stars. As a class their spectra 
display strong Balmer emission lines and they additionally exhibit significant infrared excess rela- 
tive to the stellar photospheric continuum, attributed to emission from circumstellar gas and dust 
([Hernandez et al.ll2005l ; IUemura et alj|2004l ). Herbig Ae/Be stars are further characterized by aperi- 
odic large-amplitude (up to 4 mag) variability with minima that initially re dden and subsequently 
undergo a color "turna round" as the star becomes bluer as it gets fainter (jHerbst Shevchenko 



1999 



Grinin et allboOlh . They are catagorized within the class of UXors (named afte r UX Orion) 



since it is thought that there is n o intrinsic difference between the two (|Herbstlll994! : iNatta et al 



1997 



Herbst fc Shevchenkolll999h . 



We flag 20 of our variables as potential Herbig Ae/Be stars. Eight of these show large-amplitude 
(>2.0 mag) variations while five of the remaining vary with lower amplitudes (<1.0 mag). Figure[7| 
shows a histogram of -R-band magnitude ranges for all 20 potential Herbig Ae/Be stars. The most 
interesting features are the dip in amplitude in the medium amplitude range (1.5-2.0 mag), the 
subsequent rise in amplitude (2.0-2.5 mag), and the presence of multiple (two) large-amplitude 
(>3.0 mag) variables. Considered as a whole, the amplitud e distribution for these poten tial Herbig 
Ae/Be stars in Cyg OB2 is in agreement with that found by lHerbst fc Shevchenko! (|1999l ). including 
the three aforementioned features. 



Several of our candidates display variations strikingly similar to known Herbig Ae/Be stars. 
Four follow the pattern of MiSV1147 (jUemura et al.ll2004l ). in which there are two distinct states. 
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The stars exhibit a calm state with variations of small (< 0.5 mag) amplitudes and also experi- 
ence fading episodes with deep minima of ~1 (Varid 041), ~2 (Varids 087, 100), and even ~3+ 
(Varid 116) magnitudes. Figure [8] shows a two-band light curve of such an object. Many of the 
remaining candidates mor e closely resemble RR T au ' spending roughly equal amounts of time at 
different magnitude levels (jHerbst fc Shevchenkolll999l ). Figure [9] shows an example light curve of 
a continuously active star. 

The easiest way to confirm these as Herbig Ae/Be stars would be to obtain spectra that 
display Ha emission, as well as other Balmer emission lines. Spectral classification that would 
identify these as K0 or earlier would at the least confirm them as intermediate-mass stars, though 
the lack of periodic signals that we find in their light curves is consistent with the results of 
Herbst fc Shevchenkd (119991 ) and perhaps indicates that these candidates similarly lack the surface 
hot or cool spots that are common among lower-mass PMS stars and that evince themselves as 
periodic modulations in the light curves. A population of Herbig Ae/Be stars embedded within the 
central core of Cyg OB2 would help to constrain the age of the cluster. Without further spectral 
information, however, the true status of these stars remains unconfirmed. 



5. Catalog Matching 
5.1. 2MASS Data 



Cyg OB2 suffers from high and variable extinction, Ay ~ 5-20 magnitudes ([Knodlseder 



2000 ), that is largely mitigated when observing in IR or NIR wavelengths. We extract all 2MASS 
(jSkrutskie et al.1 120061 ) point sources within 30' of the center of our survey and remove any with 
errors >0.2 mag or a null value in any band (J, H, or K), leaving 25,012 total sources. Using a 
matching radius of 1" we find unique matches to 99 (~82%) stars in our variability catalog. Fig- 
ureHUshows a CMD for all 2MASS point sources remaining after the quality cut, with the different 
types of variability labeled in different shapes and colors. 



5.2. Chandra Catalog 



The high extinction in the regio n can also be obviat ed with X-ray observatio ns, as X-rays can 
penet rate extinction up to Ay ~ 500 (jGrosso et al.ll2005l ) . The Chandra catalog of lWright &: Drake 
(J2009J) contains 1,696 X-ray sources detected within Cyg OB2, with <1% of false sources. A total 
of 97 of our variables lie within one of their two Chandra fields, and using a matching radius of 1" 
we obtain unique matches to 49 (~51%) of those 97 stars in our variability catalog. They employ a 
one-dimensional Kolmogorov-Smirnov test to investigate X-ray variability, and of our 49 matched 
stars, seven display X-ray variability with three more being possibly variable. 



They find optical and NIR counterparts for 1,501 of their sources. Using a (J, J-H) CMD 
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they calculate individual stellar masses by tracing extinction vectors. Figure [T2l shows X-ray flux 
plotted as a function of the estimated stellar mass, with stellar variability classes overplotted in 
color. It is clear that their mass estimates are reasonable for the six OB binaries, shown in purple, 
all of which have masses >15 M & . Eight of our EBs and EB-candidates that match and have mass 
estimates appear to be low-mass (<few Mq), while two are of intermediate (2-8 M®) mass. The 
sample of our EBs and EB-candidates with mass estimates is incomplete, but such a diagnostic 
would be useful for constraining the binary fraction among high-mass early-type stars. 



6. Summary 



In this work we present a catalog of photometrically variable stars in the Cygnus OB2 associ- 
ation. After g enerating R- an d /-band light curves using ISIS, we run AoV and LS period searches 
and apply the IStetsonl ([1996) variability index, finding 121 stars with two-band variability, 116 of 
them new. Of these, 99 have unique 2MASS NIR matches and 49 have Chandra X-ray matches, 
with seven exhibiting X-ray variability. We find 10 eclipsing binaries (EBs), 17 EB-candidates, and 
20 potential Herbig Ae/Be stars, with most of the rest constituting an assortment of short and 
medium period pulsating variables. 

Follow-up work is underway to obtain spectra of the potential Herbig Ae/Be stars, as Ha 
emission would be conclusive proof of their classification and help constrain the age of Cyg OB2. 
Additionally, obtaining spectral types of the EBs and EB-candidates would be useful for constrain- 
ing the binary fraction of massive stars. Cyg OB2 has the largest number of spectroscopically 
i dentified O stars (currently 65) as well as a high concentration of as many as 2,400 OB stars 
(jKnodlsederl l2000i ) . making it an excellent target for finding massive binaries. Additionally, ac- 
curate mass determinations of these potential early-type stars could add to the total number of 
known OB stars in Cyg OB2 and also provide evidence for or aga inst the proposed propensity of 
massive binaries to form as "twins", with a mass ratio near unity (jPinsonneault Stanekl 120061 ) . 
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Fig. 1. — A red filter DSS image of the Cygnus OB2 region. The image size is 2 x 2 deg 2 and 
centered on (R.A., Dec) = (20:33:12, +41:19:00) with north up and east to the left. The area 
covered by this study is shown in red whil e the variability study o f iPigulski &: Kolaczkowskil (jl998) 
is shown in green and the X-ray study of I Wright fc Drake! (|2009i ) is shown in dashed blue. 
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Fig. 2. — The rms as a function of magnitude for the 23,561 i?-band light curves (top panel) 
and 12,611 /-band light curves (bottom panel). The 121 stars that display both R- and /-band 
photometric variability are shown in red. 
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Fig. 3. — Two-band light curves of four EBs spanning a representative magnitude range. 

Table 1. Summary of Observations 



Field 


02000 


feooo 


Exp Time [s] 


No. 


of Frames 








R 


1 


R 


I 


Field 1 


20 33 16.91 


+41 14 14.5 


60 


10 


275 


228 


Field 2 


20 31 52.93 


+41 27 00.4 


60 


10 


293 


232 



-14- 



VarID = 039 Period = 2.258d 



,— , 16.8 
S 16.85 



16.9 
16.95 
15.75 
15.8 
£ 15.85 
^ 15.9 
15.95 



17.75 
17.8 

cd 

17.85 
OS 17.9 
17.95 

_ 16 " 4 
g> 16.45 

16-5 

™ 16.55 

16.6 



J 1 1 1 1 1 1 


i i i i i i i i i i |_ 

» JL • * ~ 


~~ • fit. ** 

- ™ ** 

: i i i i i i 


it V i 

i i i i i i i i i : 


* » * 

- * • • *** *« 

£ 

: • 

-I i i i i 1 i i 


* * ~ 

* * * * 

* il ill mM * 

* • _ 

1 1 1 1 1 1 1 1 1 1 T 1 F 


0.5 
VarID = 067 


1 1.5 2 
Phase 
Period=1.699d 


_i i i i i 1 i i 

_ * 


1 ! 1 1 1 1 ! 1 1 ! 1 1 j. 

\ JFK i 

» * •- 


: i i i i i i 

T i i i i 1 i i 


i i i i i i i i i : 

1 1 1 T 1 T 1 1 1 1 T 1 T 



0.5 1 1.5 
Phase 



16.9 
U 17 

a 

— 17.1 
OS 

17.2 
15.4 

'"ao 15.5 
cd 

£ 156 

~ 15.7 

15.8 



VarID = 047 Period = 0.878d 
1 — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — p 



Ay ! * Ay ! l 



t 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

* 

A' A," '• 



J I I I I I I I I I I I I I I I I I L 



0.5 1 1.5 2 

Phase 

VarID = 093 Period = 2.462d 



18.9 
Jl 19 
OS 19.1 
19.2 
17.4 



% 17.6 



„ 17.8 
18 



i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — r 



It 



I I I I I I I I I I I I I I I 



* * 



J I I I I I I I I I I I I I I I I I L 



0.5 1 1.5 
Phase 



Fig. 4. — Two-band light curves of four EB-candidates. 
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Fig. 5. — Two-band light curves of four pulsating variables. 
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Fig. 6. — A histogram of the period distribution of the 49 pulsating variables with periods <20 
days. Not shown are three remaining pulsating variables that have periods of r*22, 27, and 59 days. 
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Fig. 7. — A histogram of i?-band amplitudes for the 20 potential Herbig Ae/Be stars. 
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Fig. 8. — The light curve of a potential Herbig Ae/Be star that is similar to that of MisV1147 
(|Uemura et al.l 12004 ). exhibiting two distinct states. One is a relatively calm state with lower 
amplitude (<0.5 mag) variations while the other is an active state with deep minima that can 
reach 3+ magnitudes. 
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Fig. 9. — The light curve o f a candidate Herbig Ae/Be star that is similar to that of RR Tau 
(jHerbst &: Shevchenkdll999l ). It is continuously active and spends roughly equal amounts of time 
at all magnitude levels. 
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Fig. 10. — An (I, R-I) CMD for all 11,690 stars with two-band photometry. The different vari- 
ability classes of our 121 photometric variables are overplotted in color. Blue triangles represent 
eclipsing binaries and eclipsing binary candidates (27 in total), red squares indicate potential Her- 
big Ae/Be stars (20), purple circles mark pulsating variables (52), and light green diamonds denote 
all other variables (22). 
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Fig. 11— A near-IR CMD for the 25,012 2MASS point sources within 30' of the field center 
and with magnitude errors <0.2 mag. Our 99 variables with 2MASS matches are overplotted in 
color. Blue triangles represent eclipsing binaries and eclipsing binary candidates (26 with 2MASS 
matches), red squares indicate potential Herbig Ae/Be stars (18), purple circles mark pulsating 
variables (39), and light green diamonds denote all other variables (16). 
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Fig. 12. — X-ray flux as a function of stellar mass for the 1,417 sources from the catalog of 
Wright &: Drake! (120091 ) that have stellar mass estimates. Our 44 variables with X-ray matches 



and stellar mass estimates are overplotted in color. Blue triangles represent eclipsing binaries and 
eclipsing binary candidates (10 with matches and mass estimates), red squares indicate potential 
Herbig Ae/Be stars (6), purple circles mark pulsating variables (20), and light green diamonds 
denote all other variables (8). Six previously known OB binaries additionally have X-ray matches 
and are shown in black pentagons. 
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Table 2. Previously known photometrically variable stars in Cyg OB2 



Lit ID 


Var Type 


Sp Type 


V 
[mag] 


Reference 


Matched? 
R I 


Schulte 7 


Irr 


03If 


10.55 


1 


S 


s 


Schult 


e 18 


llT 


Bllb 


11.01 


1 


S 


s 


Schulte 27 


SB2/EW/KE 


09.5V 


12.32 


1 


s 


s 


Schulte 57 


Unk 






1 


V 


V 


V0729 


Cyg 


EB 


06.5-7 + early B? 


9.21 


2.3 


S 


S 


V2185 


Cyg 


EA EB 


09V + B9V-A0V 


12.86 


1 


s 


c 


V2186 


Cyg 


EA EB 


BOV + ? 


12.98 


1 


s 


V 


V2187 


Cyg 


Ceph 




15.39 


1 


N 


N 


V2188 


Cyg 


Be 


BlVe-B3Ve 


14.88 


1.8 


S 


N 


V2189 


Cyg 


EA EB 




16.42 


1 


N 


N 


V2190 


Cyg 


P Ceph cand 


B2Ve 


14.06 


1.8 


C 


C 


V2191 


Cyg 


EA EB 


B4Ve 


14.41 


1.8 


C 


c 


MT91 


512 


Irr 




12.73 


1 


s 


N 


MT91 


390 


Irr 


08V 


12.95 


1 


s 


c 


MT91 


436 


Unk 




16.08 


1 


N 


N 


MT91 


448 


Unk 


06V 


13.61 


1 


s 


N 


MT91 


456 


Unk 




13.53 


1 


s 


C 


MT91 


460 


Irr 




15.67 


1 


N 


V 


MT91 


503 


Irr 




15.56 


1 


c 


c 


MTE91 831 


Irr 






1 


V 


V 


MTE91 849 


EA EB 






1 


V 


V 


MTE91 876 


Unk 






1 


V 


V 


MTE91 895 


Irr/L 






1 


N 


N 


MTE91 900 


Irr/L 






1 


V 


V 


MTE91 911 


Irr/L 






1 


N 


N 


MTE91 916 


Irr/L 






1 


V 


V 


MTE91 973 


EB? 






1 


c 


V 


MTE91 974 


Unk 






1 


N 


N 


MTE91 975 


Unk 






1 


N 


N 


A36 


SB2/EA 


BOIb + B0III 


11.4 


5-7 


O 


O 



Note. — O: outside our FoV; S: saturated in our data for that band; N: within our FoV 
and unsaturated but unmatched; C: matched but photometrically constant; V: matched 
and variable 



References . - (1) Pigulski & Kolaczkowski (1998, 
Linder et all j200E); (4) iRios fc DeGioia-Eastwoodl J2004 
Oterd <2008j |; (7) Hansonl J2003I) ; (8) iKiminki et al.l J2007 



(2) iRauw et alj jl999h: (3) 
(5) IKiminki et al.1 fcOOg) ; (6) 
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Table 3. Spectroscopic binary stars in Cyg 0B2 



Lit ID 


Bin Type 


Sp Type 


V 


Reference 


Matched? 








[mag] 




r\. 


T 

1 


Schultc 1 


SB1* 


08V 


11.18 


2 


S 


N 


Schulte 2 


SB1* 


B1I 


10.64 


1 


s 


N 


Schulte 3 


SB2/EA? 


06IV? + 09III 


10.35 


3 


o 


o 


Schulte 8a 


SB2 


05.51 + 06? 


8.99 


6,7 


s 


s 


Schultc 8c 


SB1? 


05III 


10.19 


1 


s 


s 


Schulte 9 


SB2 


05? + 06-7? 


10.78 


4 


s 


s 


Schulte 11 


SB1* 


051 


10.03 


1 





o 


Schulte 15 


SB1* 


08V + B 


11.10 


2 


s 


s 


Schulte 18 


SB1* 


B1I 


11.01 


1 


S 


s 


Schultc 20 


SB1* 


09III + mid B 


11.52 


3 


s 


s 


Schulte 21 


SB1* 


BOIb 


11.42 


1 


S 


s 


Schulte 26 


SB1* 


B1III 


11.78 


1 


s 


N 


Schulte 29 


SB1? 


07V 


11.91 


1 


o 





Schulte 41 


SB1*/SB2? 


BOV 


13.49 


1 


S 


c 


Schulte 73 


SB2 


08III + 08III 


12.5 


3 


o 


o 


Schulte 74 


SB1* 


08V 


12.51 


1 


o 


o 


V1393 Cyg 


SB1* 


BOIab 


11.07 


1 


S 


s 


V1827 Cyg 


SB2 


07? + 09? 




9,10 


o 


o 


MT91 021 


SB1* 


B2II 


13.74 


1 


o 


o 


MT91 138 


SB1* 


081 


12.26 


1 


s 


s 


MT91 174 


SB1* 


B2IV 


12.55 


1 


s 


N 


MT91 196 


SB1* 


B6V 


14.81 


1 


c 


N 


MT91 202 


SB1* 


B2V 


14.40 


1 


c 


c 


MT91 234 


SB1* 


B2V 


13.25 


1 


s 


c 


MT91 238 


SB1? 


B1V 


14.91 


1 


N 


c 


MT91 241 


SB1* 


B2V 


13.41 


1 


s 


c 


MT91 252 


SB2 


B2III + B1V 


14.15 


2 


c 


c 


MT91 268 


SB1* 


B2.5V 


14.38 


1 


c 


c 


MT91 292 


SB1*/SB2? 


B2V 


13.08 


1 


s 


c 


MT91 298 


SB1? 


B3V 


14.43 


1 


c 


c 


MT91 311 


SB1* 


B2V 


13.87 


1 


c 


c 


MT91 336 


SB1* 


B3III 


14.13 


1 


o 


o 


MT91 365 


SB2? 


B1V 


13.81 


1 


s 


N 


MT91 372 


EA*/SB1 


BOV + B2?V 


14.97 


3,5 


o 


o 


MT91 403 


SB1* 


B1V 


12.94 


1 


o 


o 


MT91 428 


SB1* 


B1V 


14.06 


1 


s 


c 


MT91 448 


SB1* 


06V 


13.61 


1 


s 


N 


MT91 485 


SB1? 


08V 


12.06 


1 


s 


c 


MT91 492 


SB1* 


B1V 


14.85 


1 


N 


N 


MT91 493 


SB1* 


B5IV 


14.99 


1 


o 


O 


MT91 513 


SB1* 


B2V 


14.26 


1 


s 


c 


MT91 517 


SB1? 


B1V 


13.74 


1 


s 


N 


MT91 561 


SB1*/SB2? 


B2V 


13.73 


1 


s 


C 


MT91 573 


SB1? 


B3I 


13.87 


1 


s 


N 


MT91 605 


SB2? 


B1V 


11.78 


1 


s 


C 



Table 3 — Continued 



Lit ID 


Bin Type 


Sp Type 


V 


Reference 


Matched? 








[mag] 




R I 


MT91 620 


SB1? 


BOV 


13.89 


1 


S N 


MT91 639 


SB2? 


B2V 


14.37 


1 


C C 


MT91 692 


SB2? 


BOV 


13.61 


1 


S N 


MT91 720 


SB2 


early B + early B 


13.59 


2 


o o 


MT91 759 


SB2? 


B1V 


14.65 


1 


o o 


MT91 771 


SB2 


07V + 09V 


12.06 


2 


o o 


A45 


SB2 


B0.5V + B2V7-B3V? 




3,8 


o o 



Note. — SB 1/2 denotes a confirmed spectroscopic binary, SB 1/2* a probable binary, 
and SB1/2? a possible binary; O: outside our FoV; S: saturated in our data for that band; 
N: within our FoV and unsaturated but unmatched; C: matched but photometrically 
constant 



Reference s — (lllKiminki et al.l l l2007f); (21 iKiminki et~aH [1200811; (3llKiminki et~all 
iioosl l; (41 iNaze et al.l feOOgl); (51 IWozniak et al.1 i2004ll; (6) Romand lll969ll; (71 

onl [120031) ; Ol lStroud et al.l fcoioll ; (101 lOterd j2008bf> 



|Pe Becker et al.l ll2004ll ; (81 lHansorJ [120031 ); Ol lStroud et ahl iboioll ; (101 lOterd j2008 



Table 4. Parameters of the 121 variable stars identified in this study 



VAR ID 


C*2000 




"2000 


v en -L y jj<-< 


2MASS TP) 


R 


I 


j 


H 


K 


X- 


-ray 














[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


Matrh 


A/n ri « Vilp 

V CXI 1 CX LJ 1~ 


001 


20:31:00.18 


41 


:34:30.07 




20310017+4134300 


13.58 


13.29 


12.523 


12.089 


11.979 


o 




002 


20:33:40.09 


41 


:08:29.65 


Pul 




13.76 


13.39 








N 




003 


20:33:39.61 


41 


:10:19.26 




20333961+4110192 


13.91 


11.50 


6.728 


5.296 


4.621 


N 




004 


20:33:42.13 


41: 


: 22:22. 77 


EB 


20334212+4122227 


14.01 


13.51 


12.313 


11.905 


11.675 


O 




005 


20:34:04.04 


41 


:14:43.04 


Pul 


20340404+4114430 


14.10 


12.95 


10.660 


9.999 


9.598 


o 




006 


20:31:42.01 


41 


:21:58.29 


Pul 


20314200+4121582 


14.22 


13.98 


13.298 


12.925 


12.857 


o 




007 


20:32:19.94 


41 


:33:55.08 


Pul 


20321994+4133550 


14.52 


13.88 


12.520 


12.122 


11.917 


Y 


C 


008 


20:33:56.32 


41: 


:18:17.07 


Pul 


20335631+4118170 


14.59 


13.88 


12.345 


11.928 


11.647 


N 




009 


20:33:04.37 


41 


:24:39.96 


Pul 




14.84 


14.54 








O 




010 


20:32:34.97 


41 


:35:03.36 


EB? 


20323497+4135033 


14.84 


14.10 


12.666 


12.177 


11.941 


N 




Oil 


20:33:31.45 


41: 


: 20:57.62 


Pul 


20333144+4120576 


15.23 


14.32 


12.270 


11.424 


10.651 


N 




012 


20:34:05.78 


41 


:19:33.26 






15.31 


14.84 








N 




013 


20:33:52.24 


41 


:23:21.50 


Pul 




15.32 


14.64 








O 




014 


20:33:23.06 


41 


:10:13.86 


Pul 


20332306+4110138 


15.41 


14.45 


12.378 


11.790 


11.475 


N 




015 


20:32:09.99 


41 


:18:22.87 






15.42 


14.98 








O 




016 


20:31:23.59 


41 


: 29:48. 97 


Pul 




15.60 


14.91 








Y 


c 


017 


20:31:37.32 


41 


:24:38.21 


Pul 


20313732+4124382 


15.71 


14.91 


13.129 


12.585 


12.323 


N 




018 


20:33:39.86 


41 


:21:45.19 


Pul 


20333985+4121451 


15.83 


14.93 


12.988 


12.384 


12.017 


N 




019 


20:31:25.74 


41 


: 29:15. 77 


Pul 


20312573+4129157 


15.85 


15.40 


14.376 


13.877 


13.878 


N 




020 


20:33:08.73 


41 


:18:51.55 


Pul 




15.94 


15.07 








N 




021 


20:33:08.33 


41: 


:15:43.26 


EB 


20330833+4115432 


16.08 


15.36 


13.679 


13.095 


12.850 


Y 


c 


022 


20:32:36.58 


41 


:32:18.58 


EB 


20323658+4132185 


16.09 


15.05 


13.148 


12.533 


12.239 


Y 


V 


023 


20:32:17.79 


41 


: 20:00.85 


Pul 


20321779+4120008 


16.29 


15.20 


12.891 


11.838 


11.362 


O 




024 


20:32:30.88 


41 


:18:10.75 


Pul 


20323087+4118107 


16.30 


15.66 


14.563 


13.892 


13.738 


N 




025 


20:32:37.86 


41 


:28:52.60 


EB? 


20323786+4128525 


16.32 


15.29 


13.263 


12.361 


12.047 


Y 


c 


026 


20:32:29.87 


41 


:32:10.16 


Pul 


20322986+4132101 


16.42 


15.24 


12.986 


12.010 


11.624 


Y 


c 


027 


20:33:31.78 


41 


: 20:47.35 


Pul 


20333177+4120473 


16.45 


15.41 


13.270 


12.361 


12.038 


Y 


c 


028 


20:32:41.11 


41 


.24:12.53 




20324111+4124125 


16.46 


15.79 


14.221 


13.557 


13.355 


N 




029 


20:33:50.90 


41 


:19:03.50 


Pul 


20335089+4119034 


16.50 


14.02 


9.434 


7.713 


6.944 


N 




030 


20:32:25.91 


41: 


:35:58.92 


Pul 


20322590+4135589 


16.50 


15.41 


13.279 


12.282 


11.779 


Y 


c 


031 


20:33:47.13 


41: 


:13:16.59 




20334713+4113165 


16.51 


15.82 


14.259 


13.734 


13.573 


Y 


c 


032 


20:33:21.96 


41 


:11:40.91 


Pul 


20332196+4111409 


16.58 


15.43 


13.008 


12.204 


11.816 


N 




033 


20:33:53.10 


41: 


:11:19.40 


EB? 


20335310+4111193 


16.59 


15.36 


13.823 


13.241 


13.068 


Y 


c 


034 


20:33:48.81 


41: 


:23:08.26 


Pul 


20334880+4123082 


16.67 


15.74 


13.773 


13.203 


12.896 


O 





Table 4 — Continued 



VAR ID 


<*2000 




^2000 


Var Type 


2MASS ID 


R 


I 


J 


H 


K 


X 


■ray 














[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


Match 


Variable 


035 


20:33:15.59 


41 


:21:33.58 






16.75 


15.70 








Y 


C 


036 


20:33:22.25 


41 


:17:44.60 


Pul 


20332224+4117446 


16.77 


15.55 


13.037 


12.005 


11.605 


Y 


V 


037 


20:32:51.37 


41 


:08:40.07 


Her? 


20325137+4108400 


16.84 


15.50 


12.492 


10.988 


9.706 


N 




038 


20:33:08.60 


41 


:19:15.41 


Her? 


20330859+4119154 


16.86 


15.76 


13.168 


11.970 


10.932 


Y 


C 


039 


20:32:23.76 


41 


:29:13.03 


EB? 


20322376+4129130 


16.86 


15.85 


13.843 


13.053 


12.757 


Y 


C 


040 


20:32:30.17 


41 


:34:49.27 


EB 


20323017+4134492 


16.86 


15.85 


13.854 


13.242 


12.925 


N 




041 


20:32:23.88 


41 


:20:59.79 


Her? 


20322387+4120597 


16.93 


15.77 


13.738 


12.699 


12.028 


N 




042 


20:32:25.06 


41 


:17:21.93 


Her? 


20322506+4117219 


16.95 


15.79 


13.090 


12.018 


11.311 


O 




043 


20:32:10.88 


41 


:22:05.57 




20321087+4122055 


16.96 


15.80 


13.506 


12.597 


12.193 


N 




044 


20:33:36.95 


41 


:23:21.65 


EB? 


20333694+4123216 


16.97 


15.88 


13.691 


12.803 


12.395 


O 




045 


20:32:26.23 


41 


:29:49.51 


Pul 


20322622+4129495 


17.05 


16.06 


13.995 


13.188 


12.913 


Y 


C 


046 


20:33:18.00 


41 


:10:20.67 


Pul 


20331800+4110206 


17.06 


15.94 


13.652 


12.909 


12.608 


N 




047 


20:33:49.00 


41 


:14:09.00 


EB? 


20334899+4114089 


17.07 


15.56 


12.634 


11.737 


11.163 


N 




048 


20:33:01.78 


41 


:11:10.77 


CV 




17.14 


16.59 








Y 


V 


049 


20:32:19.92 


41 


:21:11.68 


EB? 


20321992+4121116 


17.15 


15.98 


13.644 


12.677 


12.313 


N 




050 


20:33:18.27 


41 


:14:11.49 




20331826+4114114 


17.16 


15.92 


13.323 


11.933 


10.816 


N 




051 


20:32:11.73 


41 


:26:12.21 


EB? 




17.21 


16.06 








Y 


c 


052 


20:33:58.35 


41 


:20:56.26 






17.24 


16.08 








Y 


V 


053 


20:31:53.13 


41 


:33:31.20 




20315312+4133311 


17.43 


16.53 


14.754 


13.875 


13.561 


Y 


c 


054 


20:33:08.18 


41 


:06:26.55 


Pul 


20330817+4106265 


17.44 


16.12 


13.278 


12.118 


11.699 


N 




055 


20:32:29.05 


41 


:25:37.07 


Pul 


20322904+4125370 


17.44 


16.31 


14.093 


13.185 


12.784 


Y 


c 


056 


20:32:45.02 


41 


:23:49.18 


Pul 




17.48 


16.71 








N 




057 


20:33:16.01 


41 


:17:27.26 


Pul 


20331600+4117272 


17.49 


16.28 


13.671 


12.740 


12.200 


N 




058 


20:32:08.03 


41 


:27:14.22 


EB 


20320802+4127142 


17.50 


16.53 


14.622 


13.878 


13.651 


Y 


c 


059 


20:32:25.76 


41 


:28:42.51 


Pul 


20322575+4128425 


17.52 


16.42 


13.354 


12.386 


11.871 


Y 


c 


060 


20:33:29.65 


41 


:18:45.88 


Her? 


20332964+4118458 


17.60 


16.39 


13.812 


12.300 


11.073 


N 




061 


20:33:40.81 


41 


:21:14.65 


EB? 


20334081+4121146 


17.65 


16.44 


14.076 


13.117 


12.683 


Y 


p 


062 


20:32:20.31 


41 


:23:23.81 


Pul 




17.68 


16.55 








Y 


c 


063 


20:32:25.23 


41 


:24:27.65 




20322522+4124276 


17.72 


15.23 


10.379 


8.621 


7.803 


N 




064 


20:32:07.49 


41 


:24:35.49 


EB? 


20320748+4124354 


17.73 


16.61 


14.014 


13.125 


12.873 


Y 


c 


065 


20:33:27.54 


41 


:20:15.88 


Pul 




17.78 


16.66 








Y 


c 


066 


20:33:20.16 


41 


:19:43.28 


Pul 


20332016+4119432 


17.83 


16.57 


13.942 


12.822 


12.411 


Y 


c 


067 


20:31:48.69 


41 


:21:54.97 


EB? 


20314869+4121549 


17.84 


16.49 


13.844 


12.979 


12.528 


O 




068 


20:32:01.92 


41 


:31:36.74 


Pul 




17.87 


16.82 








Y 


c 



Table 4 — Continued 



VAR ID 


<*2000 


$2000 


Var Type 


2MASS ID 


R 


I 


J 


H 


K 


X 


■ray 












[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


Match 


Variable 


069 


20:33:57.42 


41:21:48.22 


Pul 




17.95 


16.81 








O 




070 


20:32:28.42 


41:20:17.48 


Her? 


20322842+4120174 


18.00 


16.75 


14.384 


13.258 


12.475 


N 




071 


20:32:26.34 


41:24:09.54 


EB? 


20322633+4124095 


18.05 


16.83 


14.540 


13.588 


13.126 


N 




072 


20:31:22.63 


41:28:33.76 


EB 


20312263+4128337 


18.13 


16.87 


14.192 


13.368 


12.954 


N 




073 


20:32:31.81 


41:27:14.99 






18.24 


17.08 








Y 


P 


074 


20:33:41.86 


41:10:20.35 


Pul 




18.25 


16.83 








Y 


V 


075 


20:32:37.87 


41:21:19.40 




20323787+4121194 


18.25 


16.88 


14.374 


13.212 


12.805 


Y 


c 


076 


20:33:23.36 


41:19:12.63 


Pul 


20332336+4119126 


18.38 


17.22 


14.807 


13.749 


13.355 


Y 


c 


077 


20:33:04.24 


41:06:51.43 


Pul 


20330423+4106514 


18.39 


16.95 


13.945 


12.741 


12.290 


N 




078 


20:33:57.05 


41:06:54.89 


EB 


20335704+4106548 


18.40 


16.99 


14.151 


13.155 


12.713 


O 




079 


20:33:28.11 


41:21:40.26 




20332810+4121402 


18.43 


17.31 


14.876 


13.899 


13.494 


N 




080 


20:31:04.79 


41:20:26.10 


Pul 


20310478+4120261 


18.45 


17.11 


14.458 


13.393 


12.940 


O 




081 


20:33:26.34 


41:12:38.85 


Pul 


20332633+4112388 


18.46 


17.13 


14.423 


13.380 


12.914 


Y 


c 


082 


20:31:53.46 


41:23:48.71 


Her? 




18.46 


17.19 








N 




083 


20:32:53.77 


41:15:13.40 


EB? 


20325377+4115134 


18.52 


17.09 


14.122 


12.961 


12.488 


Y 


V 


084 


20:33:12.73 


41:11:33.02 




20331273+4111330 


18.55 


17.12 


13.983 


12.674 


11.989 


Y 


c 


085 


20:34:08.62 


41:06:59.75 


EB 


20340861+4106597 


18.58 


17.26 


14.366 


13.461 


13.084 


O 




086 


20:32:31.04 


41:18:55.09 


Pul 


20323104+4118550 


18.58 


17.25 


14.817 


13.713 


13.282 


N 




087 


20:31:53.57 


41:25:34.16 


Her? 


20315356+4125341 


18.66 


17.67 


14.721 


13.440 


12.545 


N 




088 


20:32:32.93 


41:18:23.14 




20323292+4118231 


18.70 


17.33 


14.617 


13.422 


12.737 


N 




089 


20:32:19.01 


41:22:24.63 




20321901+4122246 


18.76 


17.48 


15.063 


13.948 


13.629 


Y 


c 


090 


20:33:02.33 


41:23:45.52 


Pul 


20330232+4123455 


18.84 


17.56 


14.775 


13.721 


13.292 


N 




091 


20:32:23.90 


41:16:48.49 


Pul 


20322389+4116484 


18.85 


17.26 


15.002 


14.018 


13.821 


O 




092 


20:32:26.80 


41:30:38.82 


Pul 


20322680+4130388 


18.88 


17.52 


14.875 


13.748 


13.455 


Y 


c 


093 


20:32:33.60 


41:21:02.56 


EB? 


20323360+4121025 


18.94 


17.55 


14.909 


13.716 


13.314 


N 




094 


20:31:16.89 


41:30:01.23 


Her? 


20311688+4130012 


18.94 


17.89 


15.651 


14.498 


13.924 


O 




095 


20:31:51.49 


41:23:54.33 


EB 


20315149+4123543 


18.97 


17.84 


15.243 


14.371 


14.091 


N 




096 


20:34:10.98 


41:08:14.63 


Her? 


20341097+4108146 


19.11 


17.68 


14.326 


12.957 


12.181 


O 




097 


20:33:04.50 


41:24:16.25 


Her? 




19.11 


17.93 








N 




098 


20:33:09.93 


41:17:56.39 


Pul 


20330992+4117563 


19.15 


17.82 


14.827 


13.538 


12.793 


Y 


c 


099 


20:31:59.75 


41:27:23.16 


EB? 


20315974+4127231 


19.17 


17.82 


15.069 


13.971 


13.645 


Y 


c 


100 


20:33:18.89 


41:14:37.15 


Her? 


20331888+4114371 


19.23 


17.62 


14.971 


13.170 


12.009 


Y 


V 


101 


20:31:18.74 


41:34:55.14 


Pul 


20311873+4134551 


19.30 


18.08 


15.279 


14.208 


13.862 


O 




102 


20:32:31.77 


41:24:31.37 


Pul 


20323177+4124313 


19.33 


18.04 


15.184 


14.084 


13.605 


N 
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Table 4 — Continued 



VAR ID 


<*2000 


52000 


Var Type 


2MASS ID 


R 


I 


J 


H 


K 


X 


■ray 












[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


Match 


Variable 


103 


20:32:55.97 


41:22:34.33 


Pul 


20325597+4122343 


19.40 


18.11 


15.452 


14.283 


13.802 


N 




104 


20:33:11.47 


41:12:56.87 


Pul 




19.48 


18.07 








Y 


C 


105 


20:33:55.21 


41:20:01.19 


Her? 


20335520+4120011 


19.49 


17.97 


14.929 


13.659 


12.862 


N 




106 


20:32:15.63 


41:20:09.69 


Her? 


20321563+4120096 


19.55 


18.07 


15.208 


14.049 


13.761 


O 




107 


20:31:47.13 


41:26:36.18 


Pul 


20314713+4126361 


19.66 


18.33 


15.081 


14.022 


13.569 


Y 


C 


108 


20:31:30.27 


41:28:49.52 


EB? 


20313027+4128495 


19.67 


18.47 


15.866 


14.850 


14.650 


N 




109 


20:32:36.81 


41:19:44.74 


Pul 


20323680+4119447 


19.69 


18.13 


15.446 


14.233 


13.764 


Y 


C 


110 


20:33:14.35 


41:12:21.89 


Her? 


20331434+4112218 


19.71 


18.27 


15.390 


13.908 


13.076 


Y 


C 


111 


20:33:11.52 


41:07:19.73 




20331151+4107197 


19.75 


18.46 


15.707 


14.532 


14.279 


Y 


C 


112 


20:34:02.52 


41:11:15.63 


EB? 


20340252+4111156 


19.76 


18.28 


15.161 


13.884 


13.366 


O 




113 


20:33:08.69 


41:21:10.25 


Her? 


20330868+4121102 


19.89 


18.47 


15.994 


14.376 


13.409 


N 




114 


20:31:56.79 


41:21:53.64 


Pul 




19.96 


18.44 








O 




115 


20:32:22.24 


41:22:45.45 


EB 


20322224+4122454 


19.99 


18.72 


16.243 


15.127 


15.069 


N 




116 


20:31:52.41 


41:26:10.28 


Her? 


20315240+4126102 


20.00 


18.88 


16.340 


15.044 


14.408 


Y 


c 


117 


20:32:57.47 


41:04:55.19 




20325746+4104551 


20.02 


18.57 


15.295 


14.060 


13.642 


O 




118 


20:33:12.84 


41:14:40.91 


Her? 


20331283+4114409 


20.11 


18.53 


15.378 


13.799 


13.104 


N 




119 


20:33:21.85 


41:09:44.55 


Her? 


20332184+4109445 


20.14 


18.54 


15.144 


13.682 


12.979 


Y 


c 


120 


20:32:25.03 


41:25:38.52 


Her? 


20322503+4125385 


20.31 


18.87 


15.877 


14.738 


14.363 


Y 


c 


121 


20:33:40.93 


41:09:36.29 




20334092+4109362 


21.14 


19.50 


16.128 


14.904 


14.239 


Y 


p 



Note. — O: outside our FoV; N: within our FoV and unsaturated but unmatched; Y: matched; V: variable in X-rays; P: possibly variable; C: 
constant 
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